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Abstract 

This paper pursues the study carried out by the authors in Stability and 
Hopf bifurcation in the Watt governor system [14J, focusing on the codimen- 
sion one Hopf bifurcations in the centrifugal Watt governor differential sys- 
tem, as presented in Pontryagin's book Ordinary Differential Equations, pH] ■ 
Here are studied the codimension two and three Hopf bifurcations and the 
pertinent Lyapunov stability coefficients and bifurcation diagrams, illustrat- 
ing the number, types and positions of bifurcating small amplitude periodic 
orbits, are determined. As a consequence it is found a region in the space 
of parameters where an attracting periodic orbit coexists with an attracting 
equilibrium. 
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1 Introduction 



The Watt centrifugal governor is a device that automatically controls the 
speed of an engine. Dating to 1788, it can be taken as the starting point for 
the theory of automatic control (see MacFarlane and references therein). 
In this paper the system coupling the Watt-centrifugal-governor and the 
steam-engine will be called simply the Watt Governor System (WGS). See 
Section 121 for a description and illustration, in Fig. [H of this system. 

Landmarks for the study of the local stability analysis of the WGS are 
the works of Maxwell JI] and Vyshnegradskii J^l- A simplified version of 
the WGS local stability based on the work of Vyshnegradskii is presented by 
Pontryagin 13^. A local stability study generalized to a more general Watt 
governor design was carried out by Denny [T and pursued by the authors in 

Enlightening historical comments about the Watt governor local math- 
ematical stability and oscillatory analysis can be found in MacFarlane ITU] 
and Denny There, as well as in we learn that toward the mid XIX 
century, improvements in the engineering design led to less reliable opera- 
tions in the WGS, leading to fluctuations and oscillations instead of the ideal 
stable constant speed output requirement. The first mathematical analysis of 
the stability conditions and subsequent indication of the modification in the 
design to avoid the problem was carried out by Maxwell ^1] and, in a user 
friendly style likely to be better understood by engineers, by Vyshnegradskii 

From the mathematical point of view, the oscillatory, small amplitude, 
behavior in the WGS can be associated to a periodic orbit that appears from 
a Hopf bifurcation. This was established by Hassard et al. in |5j and Al- 
Humadi and Kazarinoff in Another procedure, based in the method of 
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harmonic balance, has been suggested by Denny |lj to detect large amplitude 
oscillations. 

In ^3] we characterized the surface of Hopf bifurcations in a WGS, which 
is more general than that presented by Pontryagin ^S], Al-Humadi and 
Kazarinoff PP and Denny See Theorem 14 . 1 1 and Fig. El for a review of the 
critical curve on the surface where the first Lyapunov coefficient vanishes. 

In the present paper, restricting ourselves to Pontryagin's system, we go 
deeper investigating the stability of the equilibrium along the above men- 
tioned critical curve. To this end the second Lyapunov coefficient is calcu- 
lated f Theorem 14. 4|) and it is established that it vanishes at a unique point 
(see Fig. HI and E)). The third Lyapunov coefficient is calculated at this point 
(Theorem 14. 5|) and found to be positive. The pertinent bifurcation diagrams 
are established. See Fig. El 13 and M A conclusion derived from these dia- 
grams, concerning the region — a solid "tongue" — in the space of parameters 
where an attracting periodic orbit coexists with an attracting equilibrium, is 
specifically commented in Sectional 

The extensive calculations involved in Theorems 14.41 and 14.51 have been 
corroborated with the software MATHEMATICA 5 fH] and the main steps 
have been posted in the site [T7j . 

This paper is organized as follows. In Section |21 we introduce the WGS 
and review the Pontryagin differential equations ^Hj- The stability of the 
equilibrium points is also analyzed. This section is essentially a review of 
[ini HI HI Ell- The Hopf bifurcations in the WGS differential equations are 
studied in Sections El and |31 Expressions for the second and third Lyapunov 
coefficients, which fully clarify their sign, are obtained, pushing forward the 
method found in the works of Kuznetsov jHl El- With this data, the bifur- 
cation diagrams are established. Concluding comments, synthesizing and 
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interpreting the results achieved here, are presented in Sectional 



2 The Watt centrifugal governor system 

2.1 Differential equations for the Watt governor sys- 
tem 

According to Pontryagin [12], p. 217, the differential equations of the WGS 
illustrated in Fig. [T]are 

d (p 

di/j 2 o2 • 9 . b 

—— = c \l sm(p cosif — - sm(p ip (11 

dr I m 

1 . 

where ip G (O, |) is the angle of deviation of the arms of the centrifugal 
governor from its vertical axis Si, Q E [0, oo) is the angular velocity of the 
rotation of the flywheel D, 6 is the angular velocity of 5*1, / is the length 
of the arms, m is the mass of each ball, if is a sleeve which supports the 
arms and slides along Si, T is a set of transmission gears, V is the valve that 
determines the supply of steam to the engine, r is the time, ip = dip/dr, g is 
the standard acceleration of gravity, 6 = cQ, c > is a constant transmission 
ratio, 6 > is a constant of the frictional force of the system, / is the moment 
of inertia of the flywheel, F is an equivalent torque of the load and /i > is 
a proportionality constant. The reader is referred to Pontryagin for the 
derivation of (0) from Newton's Second Law of Motion. 
After the following change in the coordinates and time 




x = ^, y = \ -^^ z = Ci -n, T = . -t, (2) 
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Figure 1: Watt-centrifugal-governor-steam-engine system. 



the differential equations can be written as 



X 



y 



dt 

/ dy 

dt 

dz 



y 



z sm X cos X — sm x — e y 



(3) 



z' = — = a (cosx — 3) 
dt ^ 

where a>0, 0</?<l and e > 0, given by 



b I 
e = — a 
m y g 



T 1 P 1 

9 I 



(4) 



are the normalized variable parameters. Thus the differential equations Q 
are in fact a three-parameter family of differential equations which can be 
rewritten as x' = /(x, /i), where 

X = (x, y, z) G (^0, ^ j X M X [0, cx)), /i = a, e) E (0, 1) x (0, oo) x (0, oo) 

(5) 

and 

/(x, /i) = (y, sinx cosx — sinx — e y,a (cosx — /?)) . (6) 



2.2 Stability analysis of the equilibrium points 

The differential equations Q have one admissible equilibrium point 

^0 = {xo, Vo, zq) = l^arccos/3, 0, . 

The Jacobian matrix of / at Pq has the form 

/ 1 



Df (Po) 



1^ 



\ 

-e 2v//3(l-/32) 



V -a^/l-(3^ 
and its characteristic polynomial is given by p{X), with 




Figure 2: Surface of critical parameters £c = 2 a 
Theorem 2.1 For a// 

e>2a 

i/ie VFCS" differential equations (0) /iO'we an asymptotically stable equilih 
point at Pq. If0<£<2a then Pq is unstable. 



The proof of this theorem can essentially be found in Pontryagin ; it 
has also been established in a more general setting in [Tij . 

The surface of critical parameters /io = (/?, a, ^c) such that Sc = e{P, a) = 

2 a 0^1'^ is illustrated in Fig. |21 In the Section E] we will analyze the stability 
of Pq as = 2 a The change in the stability at the equilibrium Pq as 
the parameters cross the critical surface produces a Hopf bifurcation in the 
WGS, whose analysis has been carried out by PP, [3] and, in a more general 
setting, by [T^ . 

From (0)), e represents the friction coefficient of the system. The case 
£ = maybe of theoretical interest due to its connection with conservative 
systems. However, as made explicit in Vyshnegradskii's Rules, friction is 
an essential ingredient to attain stability. This point is neatly presented in 
Pontryagin , of which Figure |21 is a geometric, dimensionless, synthesis. 

3 Lyapunov coefficients 

The beginning of this section is a review of the method found in ^8j, pp 
177-181, and in for the calculation of the first and second Lyapunov co- 
efficients. The calculation of the third Lyapunov coefficient has not been 
found by the authors in the current literature. The extensive calculations 
and the long expressions for these coefficients have been corroborated with 
the software MATHEMATICA 5 
Consider the differential equations 

x' = /(x,/x), (11) 

where x G and fi G M™ are respectively vectors representing phase vari- 
ables and control parameters. Assume that / is of class C°° in M" x M*". 
Suppose (fTTj) has an equilibrium point x = xq at yU = /io and, denoting the 
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variable x — xq also by x, write 



F(x) = /(x,/io) 



(12) 



as 



F(x) = ^x + ^ 5(x, x) + ^ C(x, X, x) + ^ Z}(x, X, X, x) 
"^120 ^' ^' ^' 720 ^' ^' ^' ^' 



H i^(x, X, X, X, X, X, x) + 0(1 Ixl 

5040 ^ '''''' ^ I 



where A = /x(0, fio) and 



i?,(x,y)= 5^ 



Ci(x,y,z) = 2^ 



^0 diu d^i 



A(x,y,z,u)= V 

,,,^=1 ^0 56 



A(x,y,z,u,v) = V - 

,,fc.t^=i 



€=0 



^i(x,y,z,u, v,s) = 



a6F,(0 



^ dik Oil dir dip dig 



(13) 



(14) 



€=0 



(15) 

Xj Vk Zi Ur, (16) 
Xj Vk Zi Ur Vp, (17) 
Xj yk Zl Ur Vp Sq, (18) 



Li(x,y,z,u, v,s,t) = ^ 



J dijdikdiidirdipdiqdih 



5=0 



•^j Uk Zl Ur Vp Sq tfi, 



(19) 



for i = 1, . . . ,n. 



Suppose (xo, /io) is an equilibrium point of (fTT|) where the Jacobian matrix 
A has a pair of purely imaginary eigenvalues A2,3 = ±iuJo, > 0, and admits 
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no other eigenvalue with zero real part. Let be the generalized eigenspace 
of A corresponding to A2,3. By this is meant that it is the largest subspace 
invariant by A on which the eigenvalues are A2,3. 
Let p, g G C" be vectors such that 

n 

Aq = icjQ q, A~^ p = -iuop, {p, q) = ^Piqi = 1, (20) 

1=1 

where A~^ is the transposed matrix. Any vector y E T'^ can be represented as 
y = wq+wq, where w = {p,y) E C The two dimensional center manifold can 
be parametrized by w, w, by means of an immersion of the form x = H{w, w), 
where H : C"^ ^ R" has a Taylor expansion of the form 

H{w,w) = wq + wq + ^ -^hjkW^w'' + 0{\wf), (21) 

with hjk G C" and hjk = hkj- Substituting this expression into (fTT|l we 
obtain the following differential equation 

Hy,w' + Hi,w' = F{H{w, w)), (22) 

where F is given by (fT^ . 

The complex vectors hij are to be determined so that system ()22|1 . on the 
chart w for a central manifold, writes as follows 

w' = iuow + ^ G2iw\w\'^ + G32w\w\'^ + Gi^j,w\wf + 0(|w|^), 

with Gjk G C. 

Solving for the vectors hij the system of linear equations defined by the 
coefficients of the quadratic terms of (j22p . taking into account the coefficients 
of F in the expressions (fT!?|) and (fT^ . one has 

h^y = ~A-^B{q,q), (23) 
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h2o = {2tuJoIn-A)-'B{q,q), (24) 

where /„ is the unit n x n matrix. Pursuing the calculation to cubic terms, 
from the coefficients of the terms in ()22j) follows that 

/i3o = (SiuJoIn - A)~' [3B{q, ho) + C{q, g, q)] . (25) 

From the coefficients of the terms w'^w in (|22p one obtains a singular 
system for /i2i 

{lujoln - A)h2i = C{q, q, q) + B{q, h^o) + 2B{q, h^) - Gsig, (26) 
which has a solution if and only if 

(p, C(g, g, q) + 5(g, ho) + 25(g, /m) - G^iq) = 0. 

Therefore 

^21 = (p, C(g, g, g) + 5(g, (2za;o/n - ^)"'5(g, g)) - 25(g, ^-^^(g, g))) . (27) 

The first Lyapunov coefficient li is defined by 

/i = i Re Gai. (28) 

The complex vector /i2i can be found by solving the nonsingular (n + 1)- 
dimensional system 

iuoIn-A q \ / h2i \ / C{q,q,q) + B{q,h2o) + 2B{q,hn) - G2iq 

p j \ s j \ 

(29) 

with the condition {p, h2i) = 0. 

For the sake of completeness, in Remark 13.11 we prove that the system 
(j29|) is nonsingular and that if {v, s) is a solution of (j29j) with the condition 
{p, v) = then t> is a solution of (^^1 . 
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Remark 3.1 Write = T*^ © T'^", where T'^ and T^" are invariant by A. 
It can be proved that y G T*" if and only if {p, y) = 0. Define 

a = C{q, q, q) + B{q, /i2o) + 25(g, hu) - Gsi?- 

Let (f , s) be a solution of the homogeneous equation obtained from \29j} . 
Equivalently 

{itooln - A)v + sq = 0, {p, v) = 0. (30) 

From the second equation of i3U\} . it follows that v G T'*", and thus {iujoln — 
A)v G T*". Therefore {p, {iujoln — A)v) = 0. Taking the inner product of 
p with the first equation of i3U^} one has {p, {ioJ^In — A)v + sq) = 0, which 
can be written as {p, {iuoln — A)v) + s{p,q) = 0. Since {p,q) = 1 and 
{p, {iujoln — A)v) = it follows that s = 0. Substituting s = into the first 
equation of / T^) one has {ioj^ln — A)v = 0. This implies that 

V = aq, a G C. (31) 

But = {p,v) = {p,aq) = a{p,q) = a. Substituting a = into / T^T]) it 
follows that V = 0. Therefore {v, s) = (0, 0). 
Let (f , s) be a solution of K^) . Equivalently 

{iuoQin - A)v + sq = a, {p, v) = 0. (32) 

From the second equation of it follows that v G T'^" and thus {ioJ^Ln — 
A)v G T''". Therefore {p, {iujoln — A)v) = 0. Taking the inner product of p 
with the first equation of one has {p, {iujoln — A)v + sq) = {p, a) , which 
can be written as 

{p, {iuoLn - A)v) + s{p, q) = {p, a). 

As {p,a) = 0, {p,q) = 1 and {p, {ioJoLn — A)v) = it follows that s = 0. 
Substituting s = into the first equation of results {ioJ^In — A)v = a. 
Therefore v is a solution of / f^) . 
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The procedure above will be adapted below in connection with the deter- 
mination of and h^^. 

From the coefficients of the terms w"^, w^w and w'^w'^ in (|22|) . one has 
respectively 

h^o = (4it^o4 - ^)"^ [3S(/i2o, h2o) + 4:B{q, hso) + 6C{q, g, /lao) +D{q, q, q, q)], 

(33) 

hi = {2iUoIn - A)-\W{q, hi) + B{q, ho) + 3S(/i2o, hi) 

+3C{q, q, hi) + 3C(g, q, ho) + D{q, q, q, q) - 3G2iho], (34) 

h2 = -A~^ [D{q, q, q, q) + 4C(g, q, hi) + C{q, q, ho) + C{q, q, ho) 

+2B{hu hii) + 25(g, hi) + 25(g, hi) + B{ho, ho)i (35) 

where the term —'2,hi{G2i + G2i) has been omitted in the last equation, since 
G21 + (521 = as li = 0. 
Defining 7-^32 as 

ns2 = QB{hi, hi) + B{ho, ho) + SBChi, ho) + 3B{q, h2) 
+2B{q, hi) + 6C(g, hu hi) + SC{q, ho, ho) + 3C(g, g, hi) 
+6C(g, q, hi) + 6C(g, ho, hi) + C{q, q, ho) + D{q, q, q, ho) (36) 
+QD{q, g, g, hi) + ^D{q, g, g, /i2o) + E{q, g, g, g, g) 

—QG2ihi ~ 3^21^21) 

and from the coefficients of the terms 'ur'iiP' in (j22j) . one has a singular system 
for hi 

{iujoln - A)h2 = 'H32 - Ggag, (37) 
which has solution if and only if 

(p,H32-G'32g) = 0, (38) 

12 



where the terms —6(721^21 — 3(521/^21 in the last hne of (jHUI) actually does not 
enter in last equation, since (p, h2i) = 0. 

The second Lyapunov coefficient is defined by 



where, from (jHEl), ^32 = (p, 7^32)- 

The complex vector /i32 can be found solving the nonsingular (n + 1)- 
dimensional system 



with the condition {p, /i32) = 0. 

From the coefficients of the terms w^w, w'^w'^ and w^w^ in (j221); one has 
respectively 



/i4i = (32c^o/n - A)-^[AB{hi, ho) + 6B{ho, h2i) + AB{q, hi) 

+B{q, ho) + 12C(g, hi, ho) + 6C(g, g, hi) + 4C(g, g, ho) (41) 
+3C(g, /i2o, /i2o) + 4/^(1?, q, q, hi) + 6D{q, q, q, ho) 

+E{q,q,q,q,q) - QG2iho], 

h2 = {2iuJoIn - A)-^[8B{hi, hi) + QB{ho, ^^22) + B{ho, h^o) 
+6B{hi, hi) + ^BChi, ho) + ^B{q, h2) + 2B{q, hi) 
+12C(/in, hi, ho) + 3C(/i2o, ho, h2o) + 2AC{q, hi, hi) 
+ 12C(g, ho, hi) + 4C(g, ho, ho) + 6C(g, g, h2) + 8C(g, g, hi) 

+8C{q, hi, ho) + l2C{q, ho, h2i) + C{q, g, ho) (42) 
+ 12Z:'(g, g, hi, hi) + GD{q, g, ho, ho) + 4-D(g, g, g, hi) 
+12D(g, g, g, /i2i) + 24D(g, g, /in, ^^20) + 4Z^(g, g, g, /i3o) 




32, 



(39) 




(40) 
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+3D{q, q, /i2o, ^20) + E{q, q, q, q, hio) + SE{q, q, q, q, hu) 
+6E{q, q, g, g, /iso) + K{q, q, q, q, q, q) 
-4(6-32/^20 + 3^21/131 + 621/131)], 

/i33 = -A-^[9B{hn, h22) + 3S(/i2o, ^31) + 8^(^20, /^si) + 95(/i2i, ^21) 

+S(;i30, /i30) + 35(g, /l32) + 35(g, /l32) + 6C(/lii, /in, /ill) 

+9C(/iii, /i2o, h2o) + 18C(g, /in, /?^2i) + 3C(g, /i2o, /^so) 
+9C(g, h2o, h2i) + 3C(g, g, /131) + 9C(g, g, /i22) + 18C(g, /in, /121) 
+9C(g, /120, /i2i) + 3C(g, /120, /i3o) + 3C(g, g, /131) + 9D{q, q, /120, /in) 
+D{q, q, q, /130) + 9D{q, q, q, /121) + 18D{q, q, hu, /in) 
+9i:'(g, q, /120, /i2o) + 9D(g, q, q, /121) + 9D(g, q, hu, /120) 
+3-E(g, g, q, q, ho) + 9E{q, q, q, q, hu) + ^E{q, q, q, q, /120) 
+K{q, q, q, q, q, q) - 8(632 + G'32)/iii - 9(^21 + G'2i)/i22]. 

Defining Has as 

7^43 = 12S(/lii, /I32) + 65(/l20, /132) + 35(/l20, /i4l) 

+185(/i2i, /122) + 125(/i2i, /131) + 45(/i3o, /i3i) + BChso, h^o) 
+4S(g, /133) + SB{q, /142) + 366(/iii, hu, /i2i) + 366(/iii, /120, ^21) 
+126(/iii, /120, /iso) + 36(/i2o, /i2o, /i3o) + 186(/i2o, ho, h2i) 
+36C{q, hu, h22) + 126(g, /120, /^si) + 126(g, /120, /^si) 
+366(g, /121, ^21) + 46(g, /130, ^30) + 66(g, g, ^32) 
+126(g, q, /132) + 246(g, /in, /131) + 186(g, /120, /i22) 
+SC{q, ho, /i4o) + 186(g, /121, /121) + 126(g, ^121, /130) 
+3C(g, g, /?.4i) + 2AD{q, hu, hu, hu) + 36Z^(g, hu, /i20, /i2o) 
+36D(g, q, hu, hi) + QD{q, Q, /i20, ^30) + 18-0(5, q, ho, hi) 
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+4D(g, q, q, li^i) + 18D(g, g, g, /i22) + 72D(g, g, /in, /i2i) (44) 
+36D(g, g, /i2o, /i2i) + 12£'(g, g, /i20, /^ao) + ^'^D{q, q, q, h-n) 
+36D{q, /ill, hu, /120) + 9L'(g, /i2o, ^20, ^20) + 12/^(^, q, hi, ho) 
+18D{q, g, /i2o, /i2i) + ^(g, g, g, /i4o) + 12-E(g, g, g, /in, /i2o) 
+E(g, g, g, g, /130) + 12E(g, g, g, g, /i2i) + 36E(g, g, g, /in, /in) 
+ 18E(g, g, g, /120, /i2o) + 18-E(g, g, g, g, /121) + 36E(g, g, g, hu, ho) 
+4E(g, g, g, g, /130) + 3^(g, g, g, /120, /i2o) + 3i^(g, g, g, g, g, /120) 
+ 12K{q, g, g, g, g, hu) + QK{q, g, g, g, g, /i2o) + ^(g, g, g, g, g, g, g) 

-6(26-32/^21 + G32/121 + 3G2l/l32 + 2^21/132), 

and from the coefficients of the terms w^w^, one has a singular system for 

h43 

{iuoln - A)hi^ = H43 - G43q (45) 

which has solution if and only if 

(p,7^43-G43g) = 0, (46) 

where the terms —6(2^32/121 + G32hi + 3(^21^32 + 2G21/132) appearing in the 
last line of equation (ji^ actually do not enter in the last equation, since 
(P, ^21) = and {p, /132) = 0. 

The third Lyapunov coefficient is defined by 

^3 = ^ Re G43, (47) 

where, from (jlH)), G43 = (p, 7^43). 

The expressions for the vectors ho, h^o, /151, ho, /igi, /i52 have been omit- 
ted since they are not important here. 
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Remark 3.2 Other equivalent definitions and algorithmic procedures to write 
the expressions for the Lyapunov coefficients lj,j = 1,2,3, for two dimen- 
sional systems can be found in Andronov et al. and Gasull et al. 
among others. These procedures apply also to the three dimensional systems 
of this work, if properly restricted to the center manifold. The authors found, 
however, that the method outlined above, due to Kuznetsov 13 \9^, requiring 
no explicit formal evaluation of the center manifold, is better adapted to the 
needs of this work. 

A H op f point (xo,/Uo) is an equilibrium point of ()11|) where the Jacobian 
matrix A = /x(xo,/io) has a pair of purely imaginary eigenvalues A2,3 = 
±iuJo, ujq > 0, and admits no other critical eigenvalues — i.e. located on the 
imaginary axis. At a Hopf point a two dimensional center manifold is well- 
defined, it is invariant under the flow generated by (fTT|) and can be continued 
with arbitrary high class of differentiability to nearby parameter values. In 
fact, what is well defined is the oo-jet — or infinite Taylor series — of the 
center manifold, as well as that of its continuation, any two of them having 
contact in the arbitrary high order of their differentiability class. 

A Hopf point is called transversal if the parameter dependent complex 
eigenvalues cross the imaginary axis with non-zero derivative. In a neighbor- 
hood of a transversal Hopf point — HI point, for concision — with h the 
dynamic behavior of the system ()11|). reduced to the family of parameter- 
dependent continuations of the center manifold, is orbitally topologically 
equivalent to the following complex normal form 

w' = {rj + iu)w + liw\w\'^, 

w & C, 7], uj and li are real functions having derivatives of arbitrary high 
order, which are continuations of 0, c^o and the first Lyapunov coefficient at 
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the HI point. See 0. As Zi < (/i > 0) one family of stable (unstable) 
periodic orbits can be found on this family of manifolds, shrinking to an 
equilibrium point at the HI point. 

A Hopf point of codimension ^ is a Hopf point where li vanishes. It is 
called transversal if rj = and li = have transversal intersections, where 
7] = ri{fi) is the real part of the critical eigenvalues. In a neighborhood 
of a transversal Hopf point of codimension 2 — H2 point, for concision — 
with ^2 7^ the dynamic behavior of the system (jllj) . reduced to the family 
of parameter-dependent continuations of the center manifold, is orbitally 
topologically equivalent to 

w' = {rj + iLjQ)w + Tw\w\'^ + l2w\w\'^, 

where r) and r are unfolding parameters. See 0j. The bifurcation diagrams 
for ^2 7^ can be found in [8 , p. 313, and in p[^. 

A Hopf point of codimension 3 is a Hopf point of codimension 2 where I2 
vanishes. A Hopf point of codimension 3 is called transversal if 77 = 0, /i = 
and ^2 = have transversal intersections. In a neighborhood of a transversal 
Hopf point of codimension 3 — H3 point, for concision — with h ^ the 
dynamic behavior of the system (fTT|) . reduced to the family of parameter- 
dependent continuations of the center manifold, is orbitally topologically 
equivalent to 

w' = {r] + itUo)w + Tw\w\'^ + uw\w\'^ + l3w\w\^, 

where t], t and u are unfolding parameters. The bifurcation diagram for 
I3 can be found in Takens [T^ . 

Theorem 3.3 Suppose that the system 

x' = /(x, /i), X = (x, 2;), /i = a, e) 
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has the equilibrium x = for /i = with eigenvalues 

A2,3(At) = Vil^) 

where uj{0) = ujq > 0. For fi = the following conditions hold 

V{0) = 0, /i(0) = 0, /2(0) = 0, 

where h{fi) and hif^) are the first and second Lyapunov coefficients, respec- 
tively. Assume that the following genericity conditions are satisfied 

1. 7^ 0, where /3(0) is the third Lyapunov coefficient; 

2. the map fi ivif^): h{fJ'), h{fJ')) is regular at fi = 0. 

Then, by the introduction of a complex variable, the above system reduced to 
the family of parameter- dependent continuations of the center manifold, is 
orbitally topologically equivalent to 

w' = {j] + iu!o)w + Tw\w\'^ + z/wlwl'^ + l^wlwf 

where t], t and v are unfolding parameters. 

Remark 3.4 The proof of this theorem given by Takens for C°° families of 
vector fields, using the Malgrange- Mather Preparation Theorem is also 
valid in the present case of arbitrarily high, but finite, class of differentiability, 
using the appropriate extensions of the Preparation Theorem. See Bakhtin 
and Milman [W], among others. 

4 Hopf bifurcations 

The stability of the equihbrium point Pq given in ((Tj) as = ^iP, en) = 2aP^^'^ 
is analyzed here. According to (fT^ and the subsequent expressions (I14|) . (fTHj) . 
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(fTT)|) . (fTTj) . dH)) and (Uni), for Bi to Lj, one has 

/ 1 \ 



where 



-tug 2 /3 a;o 

V -a cjo y 



'1-/32 



/5 



and referring to the expressions in equations (fTT|) and (|T^ 

F(x)-Ax=(0,F2(x),F3(x)), 



where 

^2(X) 



3 2 ^3/2 2 2(2/^2 4-7/32 3 



6/3 



a; 



-4 Wo /3 ;z + (2/32 - 1) X ^2 _^ - Vi^cc'oa;^ + -/3-^/2(l - 2/32)x3z 

8 3 



31/?2 - 16 



2-4/32 



^ -x^ + -/3a;ox''z + 
120/3 3^ " 3 



3^2 



-X 2;' 



V^cuox^ + ^/3-i/^(2/3^ - l)x\ + ^/3^/^c.oxS2 



64- 127/32 7 8 ^ , 4/3 



X /3c(Jox 2; + 

5040/3 45'^ " 15 



xV + 0(||x||«), 



1 \ 1— 1 1 ^ 

F3(x) = --a/3x^ + -aV/^ ^Q^^ + 7T7"/^^^-T7^" a//^ ^^o^;^ 

2 6 24 120 

From equations ((13, (HI), (HHI), (HEl), (HH), (HE)), (HHI) and ^ one has 



y) = (0, 52(x, y), -a /3 xi 
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where 

2 (2 — 1) 

-B2(x, y) ^ -3 ujQ ^/p xi yi + 2 uq 0^^'^ H -j= {xi ys + xs yi), 

C(x, y, z) = (o, C2(x, y,z),a^/puJo xi yi zi^ , (52) 

where 

4 - 7/5^ 

C2(x, y, z) = — - — xi yi zi - Suq (3 {x^ yi Z3 + xi ^3 zi + X3 yi zi) 

+2 (2/5^ - 1) {xi 2/3 Z3 + 2:3 yi Z3 + 2:3 ys zi), 

D{x, y, z, u) = (0, D2{x, y, z, u), a/3a;i|/i2;iMi) , (53) 

where 

D2(x, y, z, u) = IbuJof^^^'^xiyiZiUi + 8 ( ^^^^ J (^a;i2/i2;iM3 + xiyi^rsWi 

+xiy32;iiii + x^yiZiUi^ - SujqP^^'^ (xiyiz^us + Xiyg^iitg 
+Xiy3Z3Ui + X3y3ZiUi + XsyiZiUs + 2:31/1 2;3Mi) , 

E(x, y, z, u, v) = (0, E2(x, y, z, u, v), -auoP^^^xiyiZiUiVi) , (54) 

where 

^ 31/32 - 16 ^ / 

£;2(x,y,z,u, v) = xiyiziuivi + 32a;o/5 [xiyiZium 

+XiyiZiU3Vi + XiyiZ3UiVi + Xiy3ZiUiVi + XsyiZiUiVi^ 
+8(1 - 2(3^) (xiyiZiUsVs + XiyiZ3UiV3 + XiyiZ3U3Vi 

+Xiy3ZiUiV3 + Xiy3ZiU3Vi + Xiy3Z3UiVi + X3y3ZiUiVi 

+X3yiZ3UiVi + xsyiziusvi + xsyiZiUiVs^ , 
K{x, y, z, u, V, s) = (0, X2(x, y, z, u, v, s), -aPxiyiZiUiViSi) , (55) 
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where 



^2(x,y,z,u,v, 



( 



2/32-1 



) 



(■ 



XiUiZiUiViSs 



+XiyiZiUiV3Si + XiyiZiUsViSi + XiyiZsumsi + XiyaZiUiViSi 



+X3yiZiUiViSi) + 32u;o/3^/^( xi?/i2;iMif3S3 + xiyiZiU3ViS3 + xiyi^sMifiSs 



+Xiy3ZiUiViS3 + X3yiZiUiViS3 + XiyiZiU3V3Si + XiyiZ3UiV3Si 

+Xiy3ZiUiV3Si + XsyiZiUiVsSi + XiyiZ3U3ViSi + Xiy3ZiU3ViSi 



+X3yiZiU3ViSi + Xiy3Z3UiViSi + XgyiZgMif iSi + X3y3ZiUiViSi 



L(x, y, z, u, V, s, t) = (0, L2(x, y, z, u, v, s, t),uJoP^^^XiyiZiUiViSiti) , (56) 
where 

L2(x, y, z, u, V, s, t) = {64uJq - 63(3)xiyiZiUiViSiti - 128a;o^ (xiyiZiUiViSit3 



+XiyiZiUiViS3ti + xiyiZiUiV3Siti + XiyiZiU3ViSiti + XiyiZ3UiViSitx + 



a^il/s^iMi^^iSiti + X3?/i2;iMif isiti j + 32/3(/3 - cc^q) (^xi?/i2;iMif iSsts + 

XiyiZiUiV3Sit3 + XiyiZiMsf iSits + XiyiZ3UiViSit3 + Xi?/3ZiMif iSits + 
X3?/iZiMii;iSit3 + XiyiZiUiV3S3ti + 2;i?/iZiM3t;iS3ti + XiyiZ3UiViS3ti + 
Xiy3ZiMit;iS3ti + X3t/iZiMif 153^1 + Xit/iZiMgf gSi^i + Xi?/iZ3Mif gSiti + 
Xiy3ZiUiV3Siti + X3yiZiMit;3Siti + Xit/iZ3M3t;iSiti + Xiy3ZiU3ViSiti + 

X3?/iZin3t;iSiti +xi2/32;3nit;isiti + X3yiZ3UiViSiti + X3y3ZiUiViSiti]. 



The eigenvalues of A (equation (jMI)) are 



Ai = -£c = -2a/3^/2, A2 = i ujQ, A3 



(57) 



and from ()2U|) one has 




(58) 



and 




(59) 
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Theorem 4.1 Consider the three-parameter family of differential equations 
The first Lyapunov coefficient is given by 

g{f3,a) = 3 + {a^ -5)f3^ + a^(3^ (61) 

is different from zero then the three-parameter family of differential equations 
(0j has a transversal H op f point at Pq for = £{P, a) = 2 a 

// {(3,a,ec) & S U U (see Fig. 0) then the three-parameter family of 
differential equations (0) has a HI point at Pq. If {P,a,ec) G S then the 
HI point at Pq is asymptotically stable and for each e < Sc, but close to Sc, 
there exists a stable periodic orbit near the unstable equilibrium point Pq. If 
{(3,a,ec) G U then the HI point at Pq is unstable and for each e > Ec, but 
close to Sc, there exists an unstable periodic orbit near the asymptotically 
stable equilibrium point Pq. 

This theorem summarizes Proposition 3.2 and Theorems 3.5, 3.6 and 3.7 
established in Equation (j61|) gives a simple expression for the sign of the 
first Lyapunov coefficient (j6U|) . Its graph is illustrated in Fig. |2l where the 
signs of the first Lyapunov coefficient are also represented. The curve /i = 
divides the surface of critical parameters into two connected components 
denoted by S and U where Zi < and Zi > respectively. 

We have the following theorem. 

Theorem 4.2 Consider the three-parameter family of differential equations 
(0) restricted to e = Ec. The second Lyapunov coefficient is given by 
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Figure 3: Signs of the first Lyapunov coefficient. 

where 

h{p, a, = -162 - 54(-9 + 37a^)p^ - 9(-126 + 61a^ + 60a^)/?^ - 18(405 
-3212a2 + ll28a'^)p^ + (13770 - 210843a2 + 113612a^ - 5533a^)/3® 

-6(2133 - 57687a^ + 38218a^ + 5186a^)/3^° 
+ (5994 - 301275a2 + 215340a^ + 284264^6 - 16022a^)/?^2 
+2(-567 + 678780^ - 45196a^ - 379430a^ + 9347a^)/?^^ 
+a2(-25029 + 95400^ + 990831a^ + 155856a^ - 21205a^)/3^^ 
+4a^(513 - a2(163340 + 1206160^ - 16768a^))/3^^ 
-2a^(-86887 - 258835^2 + 30173a^ + 7208a^)p^^ + 2a^(-96867 
-8956^2 + 23208a^)/?2^ + a^°(33671 - 582880^ - 4880a^)/?^^ 
+160^2(1503 + 718a2)/526 - 16a^^(453 + A0a^)(3^^ + 6400^^/3=^°. 
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Proof. Define the following functions 

Ti = Re(p, E{q, q, q, q, q)), T2 = Re(p, D{q, q, q, /lao)), 
T3 = Re(p, D{q, q, q, /lao)), T4 = Re(p, D{q, q, q, hn)), 
T5 = Re(p, C{q, q, hso)), Tq = Re{p, C{q, q, hi)), = Re{p, C{q, q, /i2i)), 

Tg = Re{p, C{q, h2o, /120)), Tg = Re{p, C{q, hu, hn)), 
Tio = Re(p, C{q, h2o, hu)), Tu = Re{p, B{q, h^i)), T12 = Re{p, B{q, /i22)), 
Ti3 = Re{p, B(Ji2o, hso)), Tu = Re(p, B(h2i, /i2o)), T15 = Re{p, B{hn, /i2i)). 
From (IHUj) one has 

Re G32 = Ti + T2 + 3T3 + 6T4 + T5 + 3T6 + 6T7 + STg + 6T9 + 6T10 

+2rn + 3ri2 + Tia + ST^ + 6T15. 

The theorem follows by expanding the expressions in definition of the 
second Lyapunov coefficient It relies on extensive calculation involving 
the vector q (fSHj). the vector p (fS^ . the functions B, C, D and -E, listed 
equations (jHT)), and (jMjl . respectively, the long complex vectors 

^11, ^20, ^30, ^21, ^31 and /i22, and the above functions Ti to T15. 

The calculations in this proof, corroborated by Computer Algebra, have 
been posted in jTTj. Here, the complex vectors hu, /i2o, ^30, ^21, ^31 and /i22 
have particularly long expressions. They have listed in MATHEMATICA 5 
files. 

■ 

Theorem 4.3 For the system (0) there is unique point Q = [(3,a,ec), with 
coordinates 

(3 = 0.86828033997971281542..., a = 0.85050048430685017856... 



24 



O.B 



0.6 



0.4 



0.2 



g{/S, C£) = and h{j5, a) = (i 



1 1 1 














It 










1 1 1 







0.2 



0.4 



0.6 



1 ^ 



and 



Figure 4: Signs of the first and second Lyapunov coefficients. 



Er = 1.37624106484659953171. 



where the curves h = and I2 = on the critical surface intersect and there 
do it transver sally. 

Computer assisted Proof. The point Q is the intersection of the curves 
Zi = and /2 = on the Hopf critical surface. It is defined and obtained by 
the solution of the equations 

g{[3,a) = 0, 

given in (j61|) . and 

h{(3,a) = h{(3,a,ec) = 0, (63) 

where h{P,a,ec) is given by (jHH). The existence and uniqueness of Q with 
the above coordinates has been established numerically with the software 
MATHEMATICA 5. 
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Figure in presents a geometric synthesis interpreting the long calculations 
involved in this proof. The sign of h{j3, a) gives the sign of the second Lya- 
punov coefficient (j62|) . The graph of h{(3,a) = 0, where the signs of the 
first and second Lyapunov coefficients are also illustrated. As follows, I2 < 
on the open arc of the curve li = 0, denoted by Ci. On this arc a typical 
reference point R is depicted. Also /2 > on the open arc of the curve li = 0, 
denoted by C2. This arc contains the typical reference point, denoted by T. 
See also Fig. El 

The bifurcation diagrams of the system (jH)) at the points T and R are 
illustrated in Fig. El and [71 as a consequence of P| and |T^ . 

The main steps of the calculations that provide the numerical evidence 
for this theorem have been posted in |17j . 
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Figure 5: Signs of I2 on the curve li = 0. 
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Theorem 4.4 // {(3,a,ec) G Ci U C2 then the three-parameter family of 
differential equations (0) has a transversal Hopf point of codimension 2 at 
Pq. If {(3,a,ec) € C2 then the H2 point at Pq is unstable and the bifurcation 
diagram is drawn in Fig. O // a, Ec) G Ci then the H2 point at Pq is 
asymptotically stable and the bifurcation diagram is illustrated in Fig. ^ 

This theorem is a synthesis of the discussion in the last part in the proof 
of Theorem 14.31 




Figure 6: Bifurcation diagram of the system (jSJ at point T. 

Theorem 4.5 For the parameter values at the point Q determined in The- 
orem \4.3[ the three-parameter family of differential equations ^ has a tran- 
versal Hopf point of codimension 3 at Pq which is asymptotically unstable 
since l^iQ) > 0. The bifurcation diagram of system (0j at the point Q is 
illustrated in Figs. and{^ 
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Figure 7: Bifurcation diagram of the system © at point R. 

Computer Assisted Proof. For the point Q take five decimal round-off 
coordinates (3 = 0.86828, a = 0.85050 and €c = 1.37624. For these values of 
the parameters one has 

p = {-i/2, 0.12224 - 0.31601i, 0.54878 + 0.21228i) , 

q = (-i, 0.53237, 0.79250), 
h^^ = (-1.75030,0,0.48792), 
h2o = (-2.24198 - 0.11191i, 0.11916 - 2.38715i, 0.04434 - 1.58196z) , 
h^o = (-2.68329 + 5.27951i, -8.43202 - 4.28554i, -4.24045 - 0.86409i) , 

G21 = -2.90053i, (64) 
h2i = (1.20918 + 0.65492i, -3.24920 + 0.64374i, 1.26042 + 1.11353i) , 
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Figure 8: Bifurcation diagram of the system Q at point Q. 

= (9.27690 + 25.24802z, -53.76550 + 19.75510i, -9.11345 + 11.36572i) , 
/i3i = (-25.72175 - 5.12199i, 4.47976 - 7.87822i, 6.22842 - 15.97687i) , 

= (-15.72589, 0, 10.92671) , 

G32 = -34.93331i, (65) 
= (27.17768 + 53.16361z, -57.53733 + 3.94677i, 52.73722 + 27.89259i) , 
= (-35.5370 + 180.2333i, -195.9736 - 10.0589i, -125.3480 - 33.7428i) , 
/i42 = (-778.4924 - 466.4510^, 362.1612 + 81.2385i, 390.2364 - 503.3807i) , 
/i33 = (-536.09324,0,835.33555), 
G43 = 56.23254 - 2424.27069i. (66) 
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From ((2H1), (EHI), dHI), (El, §B and © one has 

h{Q) = 0, /2(Q) = 0, h{Q) = Re 6-43 = 0.39050. 

The calculations above have also been corroborated with 20 decimals round- 
off precision performed using the software MATHEMATICA 5 ^H] • See ^\ . 

The gradients of the functions li, given in (plj) . and I2, given in (j^^ . at 
the point Q are, respectively 

(0.80095, -0.31847), (-0.38861, -0.85118). 

The transversality condition at Q is equivalent to the non-vanishing of the 
determinant of the matrix whose columns are the above gradient vectors, 
which is evaluated gives —0.80552. The transversality condition being sat- 
isfied, the bifurcation diagrams in Figs. |H1 and follow from the work of 
Takens [IF, taking into consideration the orientation and signs established 
in Theorems 14.31 and 14.41 



5 Concluding comments 

The historical relevance of the Watt governor study as well as its importance 
for present day theoretical and technological aspects of Automatic Control 
has been widely discussed by Denny jlj and others. See also (TUl Ej . 

This paper starts reviewing the stability analysis due to Maxwell and 
Vyshnegradskii, which accounts for the characterization, in the space of pa- 
rameters, of the structural as well as Lyapunov stability of the equilibrium 
of the Watt Centrifugal Governor System, WGS. It continues with recount- 
ing the extension of the analysis to the first order, codimension one stable 
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Figure 9: Bifurcation diagram of the system Q at point Ri. 



points, happening on the complement of a curve in the critical surface where 
the eigenvalue criterium of Lyapunov holds, as studied in [5] , and by the 
authors [l4j, based on the calculation of the first Lyapunov coefficient. Here 
the bifurcation analysis at the equilibrium point of the WGS is pushed for- 
ward to the calculation of the second and third Lyapunov coefficients which 
make possible the determination of the Lyapunov as well as higher order 
structural stability at the equilibrium point. See also 1313, [HI and j2] . 

The calculations of these coefficients, being extensive, rely on Computer 
Algebra and Numerical evaluations carried out with the software MATHE- 
MATICA 5 [18j. In the site [T^ have been posted the main steps of the 
calculations in the form of notebooks for MATHEMATICA 5. 

With the analytic and numeric data provided in the analysis performed 
here, the bifurcation diagrams are established along the points of the curve 
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where the first Lyapunov coefficient vanishes. Picture |H1 and IHl provide a 
quahtative synthesis of the dynamical conclusions achieved here at the pa- 
rameter values where the WGS achieves most complex equilibrium point. A 
reformulation of these conclusions follow: 

There is a "solid tongue" where two stable regimes coexist: one is an 
equilibrium and the other is a small amplitude periodic orbit, i.e. an oscilla- 
tion. 

For parameters inside the "tongue" , this conclusion suggests, a hysteresis 
explanation for the phenomenon of "hunting" observed in the performance 
of WGS in an early stage of the research on its stability conditions. Which 
attractor represents the actual state of the system will depend on the path 
along which the parameters evolve to reach their actual values of the param- 
eters under consideration. See Denny [3| for historical comments, where he 
refers to the term "hunting" to mean an oscillation around an equilibrium 
going near but not reaching it. 

Finally, we would like to stress that although this work ultimately fo- 
cuses the specific three dimensional, three parameter system of differential 
equations given by (H)), the method of analysis and calculations explained in 
Section El can be adapted to the study of other systems with three or more 
phase variables and depending on three or more parameters. 
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